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ABSTRACT 


The  development  of  a hybrid  computer  Implementation  for  the  equations 
of  motion  for  an  elastic  cable  In  two  dimensions  Is  presented.  The 
discrete  element  approach  Is  used.  In  which  the  partial  differential 
equations  that  represent  the  dynamics  of  the  cable  system  are  reduced 
to  a set  of  2N  ordinary  differential  equations  representing  N nodes. 

A time-share  procedure  Is  developed  In  which  analog  components  are  shared 
for  different  sections  of  the  cable.  The  basic  features  of  "switching" 
from  one  section  to  the  next,  and  the  logic  that  controls  it,  are  developed 
using  a simple  second  order  system  of  two  coupled  equations.  A sample 
cable  problem  (six  foot  buoy  using  four  nodes)  is  Implemented  in  full 
on  the  hybrid  computer,  and  the  results  are  compared  with  solutions 
obtained  by  an  all-digital  technique. 
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Laboratory.  The  work  request  unit  was  N68305-75-WR-60075.  The  work 
was  carried  out  under  work  unit  number  11576600. 
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INTRODUCTION 

Large  scale  cable  dynamics  problems,  when  solved  on  digital  computers, 
usually  run  at  an  Inordinately  large  ratio  of  required  computer  seconds 
to  problem  seconds.  This  Implies  costly  solutions.  Depending  upon 
the  degree  of  sophlsltlcatlon  with  which  the  cable  dynamics  are  modeled. 

It  would  not  be  uncommon  for  all-digital  solutions  to  large-scale  cable 
dynamics  problems  to  require  several  hundred  computer  seconds  per  problem 
second.  Such  ratios  do  not  encourage  design  optimization  studies  via 
computer  simulation. 

The  purpose  of  the  study  reported  herein  Is  to  explore  the  feasibility 
of  hybrid  computation  for  the  solution  of  cable  dynamics  problems,  and 
to  show  that  a hybrid  Implementation  of  cable  dynamics  problems  avoids 
the  penalty  of  long  computation  times.  The  following  points  are  addressed 
In  comparing  the  hybrid  approach  with  the  all-digital  approach. 

a.  The  accuracies  associated  with  each  approach 

b.  The  costs  Incurred  by  each  approach,  where  these  costs  are 
determined  by 

. 1)  The  relative  sizes  of  the  programming  efforts  Involved 

2)  The  relative  speed  of  computation  as  measured  by  the 

ratio  of  computer  time  to  problem  time 

c.  The  ease  with  which  parametric  studies  for  design  optimization 
purposes  can  be  handled  by  each  approach 

d.  The  Implications  of  machine  dependence  Incurred  by  each  approach 

e.  The  ease  with  which.  In  each  approach,  essentially  the  same 
programming  will  handle  different  cable  configurations 

BRIEF  DESCRIPTION  OF  CABLE  PROBLEM 

The  treatment  considered  In  this  report  follows  that  by  Ronald 
Lewis  Webster  In  "An  Application  of  the  Finite  Element  Method  to  the 
Determination  of  Non-Linear  Static  and  Dynamic  Responses  of  Underwater 
Cable  Structures”,  Cornell  University  Thesis.  It  assumes  a uniform, 
continuous,  linear,  cable  material  with  small  strains,  where  deformation 
Is  caused  by  gravity  plus  external  loads  (generally  fluid  loading). 

The  governing  equations.  In  normal  and  tangential  coordinates,  are: 
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Tg  • y/oos$-F„ 

where, 

m ■ mass  per  unit  unstretched  length 
w * weight  per  unit  unstretched  length  in  fluid  medium 
s * unstretched  length 
T * tension  In  the  cable 
E 3 Young's  modulus 
A * cross  sectional  area  of  cable 
angle  from  the  horizontal 

Ft  n * external  tangential,  normal  load  per  unit  unstretched  length 
Vt  n « tangential,  normal  cable  velocity 

To  solve  the  equations  of  motion  of  cables,  a discrete  element 
approach  Is  usually  employed.  The  effects  of  mass.  Internal  reactions, 
and  external  loads  are  lumped  at  a finite  number  of  points  (nodes), 
and  the  equilibrium  and  continuity  equations  are  applied  to  these  points, 
yielding  a discrete  set  of  equations.  This  is.  In  effect,  modeling 
the  continuous  system  as  a set  of  discrete  masses  and  massless  springs. 

The  result  Is  a set  of  dN  coupled  differential  equations,  where  N * 
number  of  nodes,  and  d * the  dimensionality  of  the  problem  (two  or  three). 

In  the  discrete  formulation,  the  partial  differential  equations 
for  a two  dimensional  problem  reduce,  for  node  1,  to  the  following  set 
of  ordinary  differential  equations: 

m^  x1  • -F*+  slnQj  - s1n91  + fJ+  cosQi  + f|_  cos9^j 

m1  *1  " +F1+  co#®1  * F1-  C0S®1-1  + s1n®1  + FI-  s1n01-l  " W1 
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where, 

= mass  associated  with  node  i 
wi  » weight  associated  with  node  i 

ft! 

F.^  s normal  ^orce  on  node  1 due  to  upper  half-element 

u 

F"  = normal  force  on  node  i due  to  lower  half-element 

F{+  = tangential  force  on  node  i due  to  upper  half-element 

fT  = tangential  force  on  node  i due  to  lower  half-element 

L.j  = distance  between  nodes  i+1  and  ij  = [<«,♦»  - »,>’  * (*i+l  • yi)?J  " 

sin0i  = (x.+1  - x.)/L1 
cosOi  = (y.+1  - yi)/Li 


The  external  forces  TN,  FT  are  generally  due  to  fluid  loading, 
with  terms  representing  added  mass,  normal  and  tangential  drag,  and 
cable  tension. 


BRIEF  DESCRIPTION  OF  SOLUTION  PROCEDURES 

Digital 

There  are  two  basic  approaches  to  solving  the  discrete  cable  problem 
on  the  digital  computer.  In  the  explicit  scheme,  the  initial  values 
of  each  time  step  are  used  to  compute  the  values  at  the  end  of  the  step, 
assuming  the  inputs  to  the  step  are  constant  within  the  time  step. 

For  an  acceptable  accuracy,  the  time  step  must  be  very  small,  leading 
to  very  long  computation  times.  In  the  implicit  scheme,  values  to  be 
solved  for  in  each  time  step  appear  in  both  sides  of  the  equations, 
which  are  then  solved  implicitly  by  iteration.  This  procedure  again 
is  very  time  consuming. 
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Analog  and  Hybrid 

In  an  attempt  to  alleviate  the  problem  of  long  computation  times, 
analog  solutions  have  been  Investigated.  If  the  number  of  nodes  Is 
small  enough,  the  entire  problem  can  be  solved  on  an  analog  computer. 

In  general,  this  Is  not  the  case;  each  node  has  three  degrees  of  freedom 
and  two  integrators  per  degree  of  freedom  (x,  x),  that  is,  six  Integrators 
per  node.  Since  a typical  analog  computer  has  on  the  order  of  40  Integrators, 
an  all-analog  problem  is  limited  to  less  than  10  nodes  just  on  the  basis 
of  required  number  of  integrators. 

In  order  to  handle  larger  problems  on  the  analog,  one  must  go  to 
hybrid  solutions.  Because  of  the  particular  nature  of  the  cable  problem, 
in  which  the  structure  of  the  equations  for  each  node  Is  the  same  (except 
for  boundaries),  a method  of  "time-  sharing"  analog  components  Is  feasible. 
Time-sharing  is  the  procedure  by  which,  at  each  time  step,  the  same 
set  of  analog  components  is  used  for  the  equations  of  successive  sections 
("blocks")  of  the  cable.  For  each  block  the  appropriate  initial  conditions 
(IC)  are  set  on  the  analog  computer  from  the  digital  computer  by  digital 
to  analog  converters  (D/A  transfers).  At  the  end  of  each  time  step 
the  final  values  are  read  from  the  analog  computer  and  stored  in  the 
digital  computer  via  analog  to  digital  converters  (A/D  transfers).  These 
final  values  are  used  as  Initial  conditions  for  that  block  In  the  next 
time  step.  An  example  of  this  procedure  Is  given  below,  where 
it  Is  assumed  that  the  equations  for  5 nodes  will  fit  on  the  analog, 
and  that  there  are  10  nodes  In  the  problem:  l.e.,  a block  consists 

of  5 nodes  and  the  problem  is  defined  by  2 blocks. 

ty  blockj  » 1)  D/A.  Set  IC's  for  nodes  1-5,  x},  1 * 1,5 

2)  Run  analog  for  At 

3)  A/D.  Read  final  values  for  nodes  1-5,  X?,  1 * 1,5 

tj  block?  - 1)  D/A.  Set  IC's  for  nodes  6-10,  x},i«  6,10 

2)  Run  analog  forAt 

3)  A/D.  Read  final  values  for  X?,  1 * 6,10 

I I 
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t2  blockj  ■ 1)  0/A.  Set  IC's  for  nodes  1-5,  X^,  1 ■ 1,5 

2)  Run  analog  for  &t 

3)  A/D.  Read  final  values  for  nodes  1-5,  X^,  1*1,5 

t2  block2  = 1)  D/A.  Set  IC's  for  nodes  6-10,  X*.  1 « 6,10 

2)  Run  analog  for  At 

3)  A/D.  Read  final  values  for  nodes  6-10,  X^,  i = 6,10 

et  cetera. 

BRIEF  DESCRIPTION  OF  A HYBRID  PROGRAM 

A hybrid  program  consists  of  several  separate  parts: 

a.  a patch  panel  on  the  analog  computer, 

b.  a logic  panel  on  the  analog  computer, 

c.  a real  time  executive  program  on  the  digital  computer, 
and 

d.  an  analog  setup  and  checkout  program  on  the  digital 
computer. 

The  equations  to  be  solved  are  apportioned  between  the  analog  and 
the  digital  computers.  A separate  "HOI"  (Hytran  Operations  Interpreter) 
program  is  written  that  checks  on  the  integrity  of  analog  components, 
compares  computed  "static  check"  values  with  measured  "static  check" 
values,  flags  wiring  errors,  and  sets  all  coefficient  devices  and  function 
generators.  The  logic  panel  is  designed  to  send  out  timing  pulses  and 
control  signals  used  for  controlling  and  synchronizing  the  analog  and 
digital  programs.  A digital  executive  program  is  written  that  does 
the  following: 

a.  the  necessary  setup  and  Initialization  of  the  computer  states 
for  real  time, 

b.  A/D  and  D/A  Input  and  output  operations, 

c.  digital  computations,  and 

d.  time  history  storage  (on  disk  or  tape),  to  be  processed 
at  the  end  of  the  real  time  run. 

In  the  programs  In  this  report,  there  Is  no  digital  computation  or  time 
history  storage,  although  these  could  be  easily  added. 
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The  procedure  used  when  actually  making  a run  Is: 

a.  mount  the  patch  panel  and  logic  panel  on  the  analog, 

b.  run  the  HOI  setup  and  checkout  program, 

c.  load  the  real  time  digital  program,  and 

d.  make  the  run. 

A schematic  of  the  logic  panel  timer,  analog  state,  and  digital 
executive  Is  shown  below.  The  analog  state  Indicates  the  condition 
of  the  analog  Integrators: 

a.  IC  (Initial  conditions):  Integrator  outputs  are  set  to 
Initial  conditions 

b.  OP  (operate  ):  Inputs  are  being  Integrated 

c.  HOLD:  Integrator  outputs  are  held  at  their  current  value. 


Logic 

Timer 


Analog 

State 


Executive 

Program 

Initialize  real  time  system 

Read  run  Input  parameters 

Load  Initial  analog  values 

Begin  real  time  loop 

Walt  for  IC  .signal  (A)  from  logic  panel 

Transfer  D/A's 

Walt  for  HOLD  signal  (C) 

Transfer  A/D's 

Load  next  D/A  settings 

Test  for  QUIT  signal  from  logic  panel 

Continue  1nftp 


In  the  hybrid  time-share  mode,  successive  timer  cycles  are  used 
successive  blocks  of  the  cable,  so  that  with  n blocks,  there  would  be 
n timer  cycles  for  each  problem  frame  time. 

There  were  two  major  steps  to  Implementing  the  hybrid  time-share 
approach.  The  first  was  to  develop  a functional  time-share  procedure. 
The  second  was  to  develop  the  hybrid  program  Itself.  These  steps  are 
discussed  In  the  following  sections. 
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DEVELOPMENT  OF  TIME-SHARE  PROCEDURE 


There  are  two  major  aspects  to  the  time  share  approach: 

a.  procedures  and  timing  control  for  switching  between 
different  blocks,  and 

b.  reliable  circuitry  and  programming  techniques  for  A/D 
and  D/A  transfers. 

The  switching  and  timing  techniques  we**e  developed  using  track/store 
(T/S)  amplifiers  as  storage  devices  Instead  of  the  digital  computer. 

Therefore  storage  could  be  accomplished  "off-line"  (not  connected  to 
the  digital  computer). 

Two  arbitrary  coupled  linear  second  order  dynamic  systems  were 
selected  as  the  problem  on  which  to  develop  this  technique.  This  choice 
was  motivated  by  the  desire  to  use  a problem  whose  solution  could  be 
obtained  analytically.  Also,  it  was  important  to  work  with  a dynamic 
system  which  could  be  easily  Implemented  "hybrid-wise"  to  prevent  complicating 
details  from  obscuring  t*e  main  purposes.  The  exact  all-analog  solution 
was  run  concurrently  with  .he  various  time-share  attempts  in  order  to 
have  a visual  aid  (strip  chart  output)  In  determining  the  success  of 
the  procedure.  4 


The  equations  programmed  were: 

*1^1  + bl*l  + clxl  + dlx2  = ® 

(1) 

a2X2  + b2x2  + C2X2  + d2xl  * ® 

(2) 

The  first  step  in  Implementing  equations  on  the  analog  is  to  re-scale 
them  to  reflect  the  maximum  values  of  each  variable.  The  equations 
are  thus  written  In  the  form: 


The  blocks  to  be  time-shared  are  the  two  separate  second  order 
systems.  This  Is  possible  because  the  same  analog  circuitry  Is  used 
to  solve  equation  (3)  as  Is  used  to  solve  equation  (4).  A summary  of 
analog  component  schematics  and  mnemonics  Is  given  In  Appendix  A,  and 
Appendix  B contains  the  analog  wiring  diagram  for  the  all-analog  solution. 
(Figure  Bl)  and  the  wiring  diagram,  logic  diagram,  and  timing  sequence 
for  the  track/store  version  of  the  time  share  procedure  (Figures  B2, 

B3,  and  B4). 

Functionally,  the  component  time-sharing  technique  can  be  described 
with  the  aid  of  the  Figures.  Referring  to  Figure  Bl,  It  Is  recognized 
that  the  diagram  In  Block  1,  as  far  as  componentry  Is  concerned.  Is 
Identical  to  that  In  Block  2.  However,  Block  1 solves  equation  (3) 
and  Block  2 solves  equation  (4),  with  the  one  difference  between  the 
two  equations  reflected  In  different  settings  on  digitally  controlled 
attenuator  (DCA)  00  (0.1}  and  OCA  13  (0.2).  The  Implementation  shown 
In  Figure  Bl  gives  solutions  continuous  In  time  for  both  Block  1 and 
Block  2.  If  discontinuous  solutions,  are  acceptable,  components  can 
be  shared.  For  example,  again  referring  to  Figure  Bl,  suppose  that 
at  time  t„  the  continuous  solutions  In  both  Blocks  1 and  2 were  frozen 
and  the  values  of  the  outputs  of  Integrators  000,  002,  010,  and  012 
were  stored.  A continuation  of  the  Block  1 and  Block  2 solutions  can 
take  place  as  follows.  The  stored  output  of  Integrator  000  Is  used 
as  an  IC  (Initial  Condition)  on  the  sane  Integrator.  The  stored  value 
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of  Integrator  002  Is  used  as  an  IC  on  Integrator  002.  The  stored  output 
of  Integrator  012  Is  used  as  an  Input  to  summer  001  through  OCA  01  (set 
at  -0.05).  The  circuit  In  Block  1 Is  now  run  for  a short  span  of  time, 

At^  (Block  1 and  Block  2 are  disconnected).  At  the  end  of  the  span, 

Atp  the  outputs  from  Integrators  000  and  002  are  stored  In  place  of 
the  previously  stored  values.  Now  using  the  same  components  as  shown 
In  Block  1 with  OCA's  00,  01,  02,  03,  and  10  set  to  the  values  for  the 
corresponding  OCA's  In  Block  2,  equation  (4)  (Block  2)  Is  solved.  With 
the  previously  stored  output  of  integrator  010  as  an  IC  on  integrator 
000, the  previously  stored  output  of  Integrator  012  as  an  IC  on  Integrator 
002,  and  the  most  recently  stored  output  of  integrator  002  as  an  input 
to  summer  001  through  DCA  01  (set  to  the  value  of  DCA  12),  the  circuit 
In  Block  1 is  run  for  another  short  span  of  time,  Atg.  At  the  end  of 
Atg  the  outputs  of  the  Integrators  000  and  002  are  stored  in  place  of 
previously  stored  values  for  Integrators  010  and  012,  respectively. 

Thus  equation  (4)  has  been  solved  for  the  small  time  span,  Atg.  Obviously, 
equation  (3)  has  not  been  solved  during  Atg,  but  has  held  its  solution 
(obtained  at  the  end  of  At^)  constant  over  the  time  span  Atg.  The 
discontinuous  solutions  continue  with  the  solving  of  equation  (3)  over 
a time  span  at^  as  follows.  With  the  output  of  Integrator  000  stored 
at  the  end  of  Atj  used  as  an  IC  on  Integrator  000,  the  stored  output 
(at  the  end  of  At^)  of  integrator  002  used  as  IC  on  integrator  002, 
and  the  value  stored  In  Integrator  012  as  an  Input  to  summer  001  through 
DCA  01  (set  to  the  value  for  DCA  POT  01)  the  circuit  in  Block  1 is  run 
for  another  short  span  of  time,  Atj.  Again,  during  Atj,  equation  (4) 
has  not  been  solved,  but  has  held  constant  the  solution  obtained  at 
the  end  of  Atg.  The  sketch  below  sumarlzes  this  technique. 
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The  solutions  obtained  by  the  component  time-sharing  techniques 
are  not  the  same  as  the  solutions  obtained  continuously;  however,  with 
the  proper  choice  of  frame  time,  the  differences  can  be  reduced  to  an 
acceptable  level. 

A more  detailed  description  of  the  track/store  storage  method  Is 
given  In  Appendix  B. 

The  next  stage  In  the  development  of  the  component  time  sharing 
technique  was  to  replace  the  track/store  (T/S)  summers  with  dlgltal- 
to-analog  (0/A)  and  analog-to-dlgltal  (A/0)  converters.  The  circuit 
Is  shown  In  Figure  B5.  The  logic  diagram  Is  the  same  as  for  the  T/S 
method  (Figure  B3),  but  In  place  of  an  analog  push  button  starting  the 
run,  the  digital  computer  sets  a control  line  high  to  start  the  run. 

In  the  diagram  of  Figure  B5,  the  Track/Store  summers  (shown  In  Figure 
B2)  are  replaced  by  A/D  and  0/A  converters.  The  Integrators  are  Initialized 
during  the  "A"  state;  the  problem  Is  run  during  the  "B"  state;  and  the 
values  of  the  Integrators  are  stored  In  the  digital  computer  during 
the  "C"  state. 

The  results  of  the  A/D,  D/A  time  share  method  used  to  solve  the 
sample  problem  are  shown  In  Figure  1 along  with  the  continuous  solution. 


HYBRID  IMPLEMENTATION  OF  FOUR  NODE  BUOY  RELEASE 

The  final  phase  of  the  feasibility  study  was  the  Implementation 
of  a sample  problem  on  the  hybrid  computer.  The  problem  chosen  was 
the  two-dimensional  buoy  release  problem,  where  a buoy  at  the  end  of 
an  anchored  cable  Is  released  from  Its  Initial  position  at  about  45 
degrees  and  allowed  to  relax  to  the  vertical  position.  The  geometry 
for  the  four  node  problem  Is  shown  In  Figure  2.  The  equations  as  programmed 
are  given  In  Figure  3. 
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FIGURE  1 - Solution  to  the  Two  Coupled  Second  Order  Dynamics 
Systems  Problem 
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x - coordinate  of  node  1 x1 
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FIGURE  2 - Geometry  of  Four  Node  Buoy  Release  Problem 


FIGURE  3 - Equations  of  Motion  for  Four  Node  Buoy  Release  Problem 


The  following  simplifications  to  the  equations  representing  the 
dynamics  of  the  buoy-cable  system  were  made  for  the  ease  of  programming. 
These  simplifications  were  made  for  this  test  case  only  to  reduce  the 
number  of  nonlinear  analog  components  required.  They  are  not  Inherent 
limitations  of  the  hybrid  approach  and  would  not  be  made  In  a normal 
simulation. 

a.  tangential  drag  was  neglected 

b.  a constant  normal  drag  coefficient  was  used  In  place 
of  a drag  coefficient  as  a function  of  Reynolds  number 

c.  In  computing  the  normal  drag,  the  unstretched  length 

of  the  cable  was  used  rather  than  the  stretched  length. 

The  first  stage  of  the  hybrid  Implementation  was  to  program  an 
all-analog  version  of  the  problem.  This  was  done  to  develop  the  analog 
circuitry  and  analog  scaling  for  the  problem.  In  addition,  the  HOI 
analog  setup  and  checkout  program  was  written,  which  could  then  be  easily 
modified  to  give  the  HOI  program  for  the  hybrid  Implementation. 

The  second  stage  was  to  develop  the  time  share  version  and  logic 
panel  using  track/store  amplifiers  as  storage  devices,  so  that  the 
procedure  could  be  debugged  "off-line". 

The  last  stage  was  to  change  the  storage  mechanism  to  the  digital, 
which  required  replacing  the  track/store  amplifiers  with  D/A  and  A/D 
converters  and  writing  the  digital  executive  program  to  control  real 
time  and  to  perform  the  I/O  Interfacing  (D/A  and  A/D). 

One  very  Important  feature  was  discovered  during  this  last  stage. 

A problem  with  drift  and  rounding  error  existed  In  using  the  A/D,  D/A 
loop.  However,  placing  a T/S  amplifier  between  all  Integrators  and 
A/D  converters  eliminated  the  problem  of  drift  prior  to  the  A/D  conversion. 
In  addition,  the  sample-and-hold  amplifiers  In  the  A/D  converters  were 
locked  Into  the  sample  mode.  With  this  latter  precaution,  the  results 
obtained  with  A/D  and  D/A  component  time  sharing  technique  agreed  exactly 
with  the  results  obtained  with  the  T/S  component  time  sharing  technique. 
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The  problem  was  Implemented  as  follows:  node  1 Is  anchored  (thus 
has  no  differential  equat1ons)and  node  4,  (the  buoy)  Is  a boundary 
point  and  Is  left  on  the  analog  computer  to  run  continuously.  The  two 
Intermediate  nodes,  2 and  3,  are  Involved  In  the  switching.  The  timing 
Is  such  that  node  4 Is  run  In  real  time,  and  nodes  2 and  3 are  each 
run  twice  real  time  for  a duration  of  half  the  time  step,  yielding  a 
frame  time  equivalent  to  that  of  node  4.  The' overall  timing  sequence 
Is: 

a.  "A"  state  (analog  IC):  D/A  transfer  for  Initial  conditions, 

b.  "B"  state  (analog  operate):  Integration  at  real  time 
for  node  4,  twice  real  time  for  node  2 or  3,  and 

c.  "C"  state  (analog  HOLD):  A/D  transfer  of  final  values 
to  be  used  as  Initial  values  the  next  time  this  node 
Is  run. 

Figure  4 shows  time  histories  of  the  x and  y displacements  for 
nodes  2,  3,  and  4 for  the  hybrid  time-sharing  solution  to  the  buoy  release 
cable  problem.  The  dots  superimposed  on  the  hybrid  results  are  the 
results  obtained  by  using  an  all-digital  computer  program.  Appendix 
C contains  the  wiring  diagrams  and  digital  computer  programs  which  make 
up  the  hybrid  simulation,  and  a listing  of  the  HOI  program  used  to  check 
out  the  hybrid  solution. 

Appendix  C,  page  C-2,  shows  the  analog  wiring  for  the  integrators 
and  the  switching  procedure.  The  wiring  is  set  up  so  that  either  the 
T/S  or  the  A/D  storage  method  may  be  used.  Function  relays  are  used 
to  switch  the  Initial  condition  Inputs  from  T/S  amplifier  to  D/A 
converter. 

It  was  found  that  by  setting  "A"  to  .003  seconds,  "B"  to  .02  seconds, 
and  HC"  to  .002  seconds,  representing  a node  frame  time  of  .025  seconds, 
the  solutions  agreed  with  the  all-digital  solution  based  upon  a node 
frame  time  of  .01  seconds. 
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RECOMMENDATIONS 


Included  below  are  several  specific  recommendations  for  design 
and  technique  modifications  to  the  hybrid  program  which  will  result 
in  improved  time  and  accuracy  and  which  can  be  easily  implemented. 

1.  Analog-Digital  I/O  Interface. 

There  are  two  types  of  hybrid  I/O  interfacing  for  A/D  and  D/A  conversion, 
(a)  a direct  Input/Output  interface  (DIO)  and  (b)  a PORT  interface. 

The  DIO  requires  the  Central  Processing  Unit  (CPU)  to  execute  instructions 
for  each  word  transferred,  whereas  the  PORT  requires  the  CPU  only  for 
the  initial  transfer.  Thus  in  programs  where  the  frametime  is  limited 
by  the  hybrid  I/O  it  is  often  possible  to  significantly  increase  the 
throughput  by  replacing  the  DIO  with  the  PORT. 

The  present  program  uses  the  DIO  interface  because  the  PORT  software 
was  not  available  at  the  time  of  this  development.  It  is  recommended 
that  in  future  studies  the  program  be  converted  to  use  the  PORT  interface. 

2.  Change  of  Position  Variables  Stored. 

It  Is  recommended  that  the  circuitry  shown  below  be  investigated 
to  replace  storage  of  the  actual  node  coordinates  with  storage  of  the 
change  in  node  coordinates. 


17 


Making  use  of  this  circuitry  Instead  of  the  circuitry  that  was 
used  In  the  hybrid  Implementation  described  In  this  report  will  result 
In  the  elimination  of  one  0/A  transfer  per  degree  of  freedom  and  the 
elimination  of  one  gain  of  10  on  amplifiers  11,  42,  3,  and  31,  thus 
reducing  the  A/D  conversion  error  by  an  order  of  magnitude. 


3.  Addition  of  Iterative  Scheme  to  Improve  Accuracy. 

The  procedure  used  In  this  report  was  one  In  which  the  position 
of  neighboring  nodes  were  held  constant  while  the  motion  of  the  node 
of  interest  was  Integrated.  Accuracy  was  obtained  by  decreasing  the 
time  step  until  the  changes  during  Lt  were  small  enough  to  yield  the 
desired  results. 

A better  procedure  would  be  to  Include  an  Iterative  scheme,  where 
several  passes  are  made  each  time  step  to  obtain  self  consistent  results. 
The  accuracy  would  be  determined  by  adjustable  convergence  criteria. 

One  such  scheme  Is  the  following: 

Suppose  we  have  a set  of  N simultaneous  differential  equations 
of  the  form: 


li  * F1  (Xi*  Xi+1*  xi-l)  x = 1.  2*  •••»  N 


The  equations  are  "relaxed"  to  solve. 


xeo 

*1 


'T+  T 


F1(X1*  XMX)>  dt  + MT> 


where. 


iS*).  % (X<k>(T  t AT)  ♦ X,  (T))  > average  value 


x|k)  • X,(T) 


Iterations  continue  until 


x|M  (T  +AT)  - xjk"1)  (T  +M)|  £ e 
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The  value  of  ( should  be  a value  realizable  with  analog  components. 

No  changes  in  the  analog  or  logic  patch  panel  would  be  required  to  implement 
this  more  general  technique.  The  only  changes  required  would  be  some 
additional  coding  in  the  digital  program. 

DISCUSSION  OF  VERSATILITY 

An  expansion  from  2 to  3 dimensions  would  be  very  straight  forward. 

The  procedures  would  be  the  same,  but  there  would  be  more  equations 

(z.  as  well  as  x.  and  y.)  and  more  terms  in  the  existing  equations. 

1 1 i*  2 2 2lk 

For  example,  L1+  becomes  [_(x1+1  -.Xj)  + (y1+1  - y^*  + Ui+1  - 

It  is  possible  that  fewer  nodes  will  fit  in  a block  and  more  blocks  will 
be  needed  per  time  step,  thus  requiring  more  computer  time  to  solve 
the  problem,  (not  as  in  the  all-digital  solution).  The  concepts,  however, 
remain  the  same  and  no  new  techniques  would  be  required. 

Parametric  studies  can  be  done  fairly  easily  with  hybrid  programs. 

The  program  is  designed  such  that  changing  parameters  from  run  to  run 
can  be  easily  accomplished  from  an  input  terminal  either  directly  or 
from  a file. 

Design  changes  would  require  more  effort.  The  initial  work  in 
general  structure,  timing,  I/O  procedures,  etc.,  would  remain  the  same, 
but  any  reprogramning  on  the  analog  would  require  development  time, 
dictated  by  the  nature  of  the  change.  The  analog  could  be  programmed 
to  handle  more  general  problems  initially,  but  at  the  sacrifice  of  analog 
components. 

Procedures  Involved  in  expanding  to  more  complex  problems  will 
depend  heavily  on  the  particular  configuration.  In  general,  more  equation 
programming  would  be  Involved,  either  on  analog  or  digital  computers. 

For  the  sample  four-node  buoy  problem,  all  of  the  calculation  was  done 
on  the  analog,  and  the  digital  was  used  solely  for  "storage"  of  current 
values  of  variables.  If  a problem  is  very  large.  It  might  be  advantageous 
to  do  some  of  the  calculations  on  the  digital  computer.  The  consideration 
then  is  whether  the  computation  time  Is  least  with  more  of  the  computation 
on  the  analog  board  (thus  fewer  nodes  per  block  and  more  blocks  per 
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frame  time)  or  with  more  computation  on  the  digital  (perhaps  requiring 

more  Iterations  for  convergence).  For  example,  one  extension  would  be 

Including  nodes  with  more  than  2 elements  attached  (e.g.,  nets).  This 

would  Involve  reprogramming  large  portions  of  the  analog  panel.  It 

might  be  decided  that  the  overall  efficiency  would  be  increased  If 

the  sines  and  cosines,  for  Instance,  where  done  digitally.  These  considerations 

would  have  to  be  worked  out  for  each  case. 

A hybrid  program  Is  fairly  machine  dependent,  In  that  every  hybrid 
facility  has  different  hardware  and  software.  The  general  capabilities, 
however,  should  be  similar  enough  that  only  a relatively  small  effort 
would  be  required  to  convert  the  hybrid  program  to  run  at  a similar 
hybrid  facility. 

CONCLUSIONS 

Both  the  track/store  time-share  and  the  hybrid  time-share  techniques 
were  tested  on  a cable  system  represented  by  four  nodes.  Results  obtained 
with  both  techniques  for  the  four  node  cable  system  were  found  to  agree 
with  one  another  and  with  results  obtained  from  an  all-analog  as  well 
as  an  all-digital  Implementation.  Although  the  test  case  was  small, 
it  Is  concluded  that  large  sophisticated  systems  can  be  successfully 
and  economically  simulated  In  real  time  via  hybrid  techniques. 

A reasonable  estimate  of  the  running  times  can  be  predicted  from 
the  test  problem.  An  Integration  Interval  of  .025  seconds  problem  time 
was  found  to  be  sufficient.  If  the  PORT  Interface  Is  used,  the  timing  for 
one  node  for  this  Integration  Interval  Is  estimated  as  follows: 

a.  200  x 10*6  seconds  are  required  for  Integrators  to  settle 
and  for  D/A  PORT  transfer, 

b.  100  x 10”6  seconds  are  required  for  operate  time,  and 

c.  200  x 10’6  seconds  are  required  for  A/D  PORT  transfer  time. 

Thus  for  each  node,  500  x 10"6  seconds  would  be  required  to  simulate 
.025  problem  seconds.  The  number  of  nodes  that  can  be  done  In  real 
time  Is  therefore  .025/(500  x 10”6)  ■ 50  nodes. 


Alternately,  It  would  be  possible  to  simulate  a cable  with  25  nodes 
In  twice  real  time,  or  100  nodes  In  *j  real  time,  or  200  nodes  In  •»  real 
time,  etc.  Thus,  on  the  basis  of  the  test  program,  it  appears  that 
a 100  node  cable  system  can  be  simulated  In  real  time  on  the  DTNSRDC 
Hybrid  computer  for  sophisticated  dynamic  situations. 

Assuming  that  the  above  conditions  are  realizable,  it  is  anticipated 
that  hybrid  computing  technology  can  have  a significant  impact  on  the 
dynamic  simulation  of  complex  cable  systems.  It  is  safe  to  predict 
that  proper  exploitation  of  this  powerful  technology  can  result  in  the 
transfer  of  dynamic  cable  system  developments  from  hardware  (expensive) 
areas  to  the  simulation  (inexpensive)  areas. 

Finally,  in  answer  to  the  points  addressed  in  the  Introduction 
concerning  the  relative  merits  of  hybrid  simulation  versus  an  all-digital 
simulation,  the  following  conclusions  are  drawn. 

a.  The  program  development  costs  for  hybrid  simulation  are 
greater  than  for  an  all-digital  simulation. 

b.  Hybrid  simulation  becomes  much  more  economical  than  all- 
digital  simulation  in  performing  production  runs.  Since 
the  hybrid  simulation  runs  real  time  or  faster,  the 
user  can  afford  to  simulate  many  cases  which  would  be 
prohibitive  in  cost  for  an  all-digital  simulation. 

c.  Different  cable  systems  may  require  separate  hybrid  programs, 
but  the  program  development  costs  would  be  small  since 

the  solution  technique  has  already  been  developed. 

There  is  little  question  in  the  minds  of  the  authors  that  a many 
node  cable  dynamics  simulation  can  be  implemented  using  the  component 
time  sharing  technique  together  with  an  iterative  scheme.  However, 
it  is  strongly  recommended  that  a moderately  sized  (20  node)  rather 
sophisticated  dynamic  cable  system  be  implemented  first;  e.g.,  a cable 
system  towed  by  a surface  ship  operating  in  a seaway.  A demonstration 
of  such  a successful  simulation  would  remove  any  uncertainties  for 
continuing  on  to  the  most  sophisticated  cable  systems  envisioned  for  the 
future. 
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APPENDIX  B 

DETAILED  DESCRIPTION  OF  TIME-SHARE  PROCEDURE 


Track  and  Store  Component  Time-Sharing  Technique 

[\]  j 

Figure®  shows  the  component  time-sharing  technique  using  T/S 
amplifiers.  The  "1"  and  "2"  states  of  amplifiers  241,  243,  and  251 
Indicate  that  the  outputs  of  those  amplifiers  are  the  Inverses  of  the 
state  "1"  Inputs  or  state  "2"  Inputs.  The  logic  variable  "first"  (see 
Figure®)  Is  used  for  putting  feedbacks  around  the  T/S  amplifiers  201, 

211,  221,  and  231.  Without  these  feedbacks,  the  T/S  summers 
will  drift  causing  their  outputs  to  be  initially  non-zero.  The 
variable  "first"  Is  also  used  to  Initialize  Integrator  202  during  the 
first  two  cycles  of  operation  since  the  T/S  amplifiers  are  Initially  zero. 

Figured  shows  a timing  diagram  for  the  logic  variables  of  Figure 
B3.  That  portion  of  the  logic  diagram  referring  to  control  line  0 
Is  also  used  for  the  A/0  and  0/A  component  time-sharing  technique 
to  be  discussed  later.  When  Push  Button  Flip  Flop  0 Is  depressed  the 
circuit  of  Figure  1 begins  to  operate.  Depressing  this  flip  flop 
also  starts  TIMER  0 (which  controls  the  component  time-sharing  techniques) 
cycling  between  "A",  "B",  and  "C".  Flip  Flop  31  Is  initially  In  the 
"1"  state,  consequently  Integrator  200  Is  Initialized  to  zero  via  T/S 
summer  201.  Integrator  002  Is  Initialized  via  DCA  111.  When  "B"  goes 
high.  Integrators  200  and  202  start  Integrating  the  equation  for  state 
"1"  (l.e.  Equation  3).  While  these  integrators  are  In  operate,  "Tj", 
remain  high  and  T/S  summers  201  and  221  track  the  output  of  Integrators 
200  and  202  respectively. 

At  the  end  of  "B",  "C"  goes  high  forcing  T/S  summers  201  and  221 
to  store  the  value  that  integrators  200  and  202  respectively  had  at 
the  end  of  "B".  While  "C"  remains  high  the  values  of  Integrators  200 
and  202  remain  at  the  values  they  had  at  the  end  of  "B".  At  the  end 
of  "C",  "A"  again  goes  high  causing  Flip  Flop  31  to  reset  making  "2" 
high.  Integrator  200  is  Initialized  zero  via  T/S  summer  211,  and  Integrator 
202  Is  Initialized  via  DCA  111.  When  "B"  goes  high.  Integrators  200 
and  202  start  Integrating  the  equation  for  "2"  (1.  e..  Equation  4). 

While  these  Integrators  are  on,  "T2"  Is  high  and  T/S  summers  211  and 
231  track  the  output  of  integrators  200  and  202  respectively.  At  the 
end  of  "B",  "C"  goes  high  forcing  T/S  summers  211  and  231  to  store 
the  values  that  Integrators  200  and  202  respectively  had  at  the  end 
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of  "B".  While  NCN  remains  high.  Integrators  200  and  202  hold  the  values 
they  had  at  the  end  of  "B".  When  "C"  goes  low,  one  complete  cycle  has 
been  completed.  For  all  remaining  cycles,  the  variable  "First"  Is  low, 
so  that  the  feedbacks  around  the  T/S  summers  are  no  longer  applied; 
and  the  Initial  condition  of  Integrator  202  will  be  determined  by  the 
state  of  amplifier  243.  This  cycle  Is  repeated  over  and  over  again  as 
long  as  Push  Button  Flip-Flop  0 remains  depressed.  In  the  mean  time  the 
circuit  of  Figure  B1  continues  to  operate. 

If  now  output  of  amplifier  002  (Figure  Bl)  and  T/S  221  (Figure  B2) 
differenced  and  the  output  of  amplifier  012  and  T/S  231  are  differenced, 
error  terms  are  generated  representing  the  difference  between  the 
continuous  solutions  and  the  solutions  obtained  by  the  time  sharing 
technique.  By  varying  the  "A",  "B",  and  "C"  time  intervals.  It  Is  possible 
to  minimize  these  errors. 


FIGURE  B1  - Analog  Wiring  Diagram  for  Continuous  Analog  Solution 
of  Two  Coupled  Second  Order  Dynamic  Systems 


Analog  Hiring  Diagram  for  Component-Tlme-Sharlng  Technique 
(Employing  Track/Store  Summers  for  Storage)  for  Solutions 
to  Two  Coupled  Second  Order  Dynamic  Systems 


FIGURE  B3  - Logic  Circuit  for  Operational  Control  of  the  Analog  Circuits 
Shown  In  Figures  1,  2,  and  5. 
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